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1 . Introduction or for correction in the visible a large number of guide stars and 

wavefront sensors are needed. 



Abstract. Laser guide stars with adaptive optics allow astronomical image correction in the absence of a natural guide star. 
Single guide star systems with a star created in the earth's sodium layer can be used to correct the wavefront in the near infrared 
spectral regime for 8-m class telescopes. For possible future telescopes of larger sizes, or for correction at shorter wavelengths, 
the use of a single guide star is ultimately limited by focal anisoplanatism that arises from the finite height of the guide star. 
^ • To overcome this limitation we propose to overlap coherently pulsed laser beams that are expanded over the full aperture of 

' the telescope, traveling upwards along the same path which light from the astronomical object travels downwards. Imaging the 

^\ , scattered light from the resultant interference pattern with a camera gated to a certain height above the telescope, and using 

phase shifting interferometry we have found a method to retrieve the local wavefront gradients. By sensing the backscattered 
light from two different heights, one can fully remove the cone effect, which can otherwise be a serious handicap to the use of 
laser guide stars at shorter wavelengths or on larger telescopes. Using two laser beams multiconjugate correction is possible, 
resulting in larger corrected fields. With a proper choice of laser, wavefront correction could be expanded to the visible regime 
and, due to the lack of a cone effect, the method is applicable to any size of telescope. Finally the position of the laser spot could 
be imaged from the side of the main telescope against a bright background star to retrieve tip-tilt information, which would 
CIh' greatly improve the sky coverage of the system. 

2 ' Key words, telescopes - instrumentation: adaptive optics - techniques: high angular resolution 

> 

■ Adaptive optics correction of wavefront distortions is becom- Another approach is to use the upward path of a projected 
' ing a common tool at large ground based telescopes. Diffraction laser beam to sample the atmospheric turbulence, either in an 
limited imaging through the turbulent atmosphere can be re- inverse Shack-Hartmann geometry by distributing many laser 
garded as routine at near infrared wavelengths and shows great beams in parallel over the telescope, or with a single laser beam 
success in many scientific areas (see for example the pro- expanded over the whole aperture of the telescope. Buscher et 
ceedings of the 'Science with adaptive optics' meeting, eds. al. I2002I have proposed a method to sample the local curva- 
Brandner & Kasper 2005 1. By using laser guide stars this tech- ture of the atmospheric wavefront distortions from the intensity 
nique can be applied at almost any position in the sky where a changes that a collimated laser beam shows when it has trav- 
suitable star for tip-tilt correction can be found. Nevertheless, eled to a certain height. In the same publication a method is 
the quality of correction when using a single guide star is lim- proposed to compare the full aperture expanded beam with lots 
ited by focal anisoplanatism, or the 'cone effect'. Due to the of other laser beams in a point-diffraction scheme. Also using 
limited height where the artificial guide star is created, some a full aperture expanded laser, Bonaccini et al. 2004 have pro- 
part of the wavefront distortion remains unsampled, and thus posed to measure the wavefront on the downward path when 
uncorrected. For possible future extremely large telescopes, or imaging back the distributed beacon with a shearing interfer- 
for correction at short wavelengths, this limit will ultimately ometer. Feasibility studies of all the above mentioned meth- 
prevent diffraction limited performance. To overcome this, to- ods are currently being carried out. Limitations for the inverse 
mographic solutions have been proposed (e.g. by Angel & Shack-Hartmann and the point-diffraction scheme may arise 
Lloyd-Hart I2000> . utilizing multiple Rayleigh guide stars or due to both the size of the individual spots that can be reached 
multiple stars in the earth's sodium layer. For a large telescope in the presence of diffraction and also beam quality aspects. For 
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the shearing measurement, the coherence length of the scat- 
tered light might be a serious issue. 

The method proposed here uses the coherent superposition 
of tilted laser wavefronts over the whole aperture of the tele- 
scope together with the method of laser phase shifting inter- 
ferometry (LPSI) in order to retrieve a local phase difference 
that can be used in an adaptive optics system in a similar way 
to the gradients retrieved from a standard Hartmann sensor. To 
adapt to different conditions, the sensing resolution can be cho- 
sen freely by adding an appropriate tilt to the beams. Detecting 
the light scattered from low altitudes will result in a ground 
layer correction. By extending the system to multiple detec- 
tions at different heights (or two beam systems) with more de- 
formable mirrors, a complete sampling of the turbulent layers 
in the atmosphere is possible. Therefore high Strehl ratios can 
be reached and multiconjugate adaptive optics correction over 
an extended field of view is possible. Additional means can 
help to expand the sky coverage further, since in a differential 
measurement far off-axis guide stars can serve as the tip-tilt 
reference. 

2. On-sky phase shifting interferometry 

Widely used for the measurement of optical surfaces, phase 
shifting interferometry (PSI) allows a phase difference to be 
measured with high accuracy. Standard laboratory interferom- 
eters make use of the comparison between a reference phase 
and a phase under test, e.g. in a Twyman-Green setup. The 
reference beam is phase shifted in several steps, and an in- 
terferogram is recorded at each step. Depending on the num- 
ber of phase steps, several algorithms are possible to calcu- 
late the phase difference and these can be found in textbooks 
(e.g. Malacara [1992> . For atmospheric turbulence probing the 
usual reference flat is of course not present, which means that 
a global wavefront can not be measured directly. 

Instead, with the scheme shown in Fig.[Q a local phase dif- 
ference can be measured that allows an interpretation as local 
gradients, and which can be treated as such for the wavefront 
correction. Assuming one pulse with two flat coherent laser 
wavefronts leaving the telescope — one slightly tilted with re- 
spect to the other — an interference pattern can be created at 
any distance. When reaching a certain height in the atmosphere, 
both waves will have accumulated local phase changes. Due to 
the tilt in the wavefront, the path through the atmosphere to the 
point (x,y,H) of each wave is slightly different. The lateral shear 
of the waves at ground level is denoted by a in Fig.[Q This cre- 
ates a phase difference <t> a which modulates the intensity of the 
interferogram that is written in the sky. With a camera gated 
to twice the time of flight to the point (x,y,H) the scattered 
light from this pattern can be imaged using the full telescope 
aperture. It should be noted that the large wavefront distortions 
present in the atmosphere result only in moderate gradients to 
be probed, and choosing the distance a smaller than the coher- 
ence length ro will avoid phase unwrapping. 

When two beams are overlayed, a tilt can only be applied in 
one coordinate direction, and a phase difference — gradient — 
can be measured in this direction. To retrieve gradients in two 
coordinates several methods are possible: such as overlapping 
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Fig. 1. Basic principle of measuring atmospheric turbulence 
with phase shifting interferometry. A pulsed laser beam is split 
into two beams which are tilted with respect to each other. To 
one of the beams a known phase shift a can be applied. The 
scattered fringe pattern which is formed in the atmosphere at 
height H, can be viewed with a camera over the full telescope 
aperture, which is gated according to the time of flight. In a 
series of laser shots a known phase difference is added to the 
beam in each pulse. Depending on the algorithm used for the 
phase retrieval, a set of three or more laser shots will allow the 
local phase difference to be measured. The frames shown here 
of a single laser shot are taken from a simulation with multiple 
turbulence layers of the fringe formation and imaging at an 8-m 
telescope with ro = 10 cm. 



more than two beams, or overlapping perpendicular polariza- 
tion measurements and disentangling the polarization planes at 
the detector. These possibilities will be outlined in Section l231 
First we describe the basic principle of a two beam overlap and 
detection. 

2.1. Formation of the interferogram 

Due to the coherence of the laser wavefronts, the intensity at 
each position (x,y) follows the equation of a two-beam inter- 
ference: 

l{x,y,d) = I'(x,y) + I"(x,y) cos(O a (x,y) + ® t (x,y) + a) (1) 

were I'(x,y) is the sum of the two beam intensities or the av- 
erage intensity and I"(x,y) represents the fringe modulation. 
<f> a (x, y) is the phase difference between the optical path lengths 
that is accumulated in the atmosphere, and <& t (x, y) denotes the 
phase difference due to the relative tilt of the two waves and all 
components present before traveling through the medium to be 
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Fig. 2. Possible arrangement for the beam combiner and phase 
shifter. The incoming laser pulse is split into two beams. In 
one path a — preferably electro-optic — phase shifter is in- 
stalled. On combining the two beams, the separation and tilt 
can be adjusted. A master clock would control the phase shift- 
ing according to the pulse sequence which is coming in. For 
high accuracy it might be necessary to have a fast control loop 
stabilizing the phase shift at the commanded value. 

probed, a is the global additional phase difference which is ar- 
tificially added to the launched beam. Shifting a in steps from 
pulse to pulse allows Q> a (x,y) to be extracted, as will be out- 
lined in Section 1231 Referring to Fig.^ it is obvious that the 
total number of fringes n resulting from the relative tilt, over 
a telescope aperture of size D, depends on the angle between 
the beams expressed here with the scattering height H and the 
chosen distance a: 



aD 
~HA~ 



(2) 



where A is the wavelength of the laser. The phase difference 
<& a (x,y) is the measurement signal that has to be extracted, and 
represents the phase gradients that have to be corrected with 
the adaptive optics system. The measured amplitude of <t> fl will 
depend on the choice of the initial shear of the laser beams and 
the relative tilt which is applied. 

When the overlapped laser wavefronts traveling upwards 
are contained in a short duration pulse, a camera which is gated 
to the time-of-flight for a certain height can image the pat- 
tern created there due to scattering by molecules. Using several 
consecutive pulses in a very short time and shifting the rela- 
tive phase of the wavefronts by a certain angle a, the unknown 
phase difference (S) a can be computed with the usual PSI algo- 
rithms. When all the intensities in Eq.^are treated as unknown, 
there are three unknown terms, and thus a minimum of three 
laser pulses are necessary to retrieve <t> fl . These pulses have to 
probe the atmosphere within a timescale in which the phase 
<t> fl can be regarded as constant. In the visible, the coherence 
time at typical astronomical sites is of order 5 ms, thus several 
pulses would need to be launched and the backscatter detected 
within this time. For a scattering height of 10 km the traveling 
time up and down is only 67 fis. But the next pulse can only be 
launched when there is no contamination from the first pulse 
on the second detection, which is approximately the travel time 
to a height of 30 km, where the last scattering occurs in the 
atmosphere, and back. This is a time of 0.2 ms, short enough 
that, including some contingency, three to four pulses could be 
launched within a millisecond. For a realistic system, a design 
optimized for the I-band and longer wavelengths would relax 
the temporal requirements, and a kHz laser and detection sys- 
tem should be fully sufficient. 



2.2. Imaging of the pattern 

The backwards imaging of the pattern with a camera will put 
some constrains on the maximum number of fringes that should 
be created. Sampling the intensity distribution, which is de- 
scribed in Eq. [2 with image points of a certain size will mod- 
ify the modulation of the detected interferogram. If the phase 
change over the pixel is linear, or nearly so, it can be shown 
easily that integration of Eq.[0over a range -A/2, A/2 leads to 
the following expression: 

I(x,y, a) — I'(x,y) 

+ I"(x,y)sinc(—)cos(® a (x,y) + ® t (x,y) + a) (3) 
2n 

with sinc(J3) = sm j^ and where A is the phase difference over 
one detection pixel. With the sine function approaching zero, 
the modulation in the interferogram approaches zero as well 
and the signal to noise is reduced. It will be shown later that 
up to A ~ 7r/2 the signal to noise is only marginally influenced. 
The main fringes are present due to the relative tilt of the wave- 
fronts. This means that with a detection camera having N pixels 
over the telescope diameter the number of fringes should not be 
larger than N/4, and therefore: 



N = 4- 



^aD 
~HA 



(4) 



To evaluate the approximate number of pixels needed we 
can choose the distance a to match one pixel, which in general 
is not a requirement. The measured <£> fl then becomes directly 
the physical phase difference from one pixel to the next and the 
number of pixels depend on the square root of 1 /HA: 



N = 



2D 

tJha 



(5) 



For an 8-m telescope, 532 nm laser wavelength and 10 km gat- 
ing this would require N = 219 elements; and for 30 km we 
find N = 126. In choosing the actual number of pixels, the re- 
duction in modulation due to the sine function has to be traded 
against the signal to noise due to the photon number incident 
on one detection element. 

On the way back from the scattering point, the light will 
be influenced by the phase distortions above the telescope, in a 
cone that has its base at the pupil and its tip at position (x,y). 
To retrieve the information written in the sky, the image of each 
point must not be distorted too strongly. A simple estimation of 
the angular distance of one element in the sky, as seen from the 
camera on the ground, shows that in general this will not be 
a serious problem. The angle is computed from H and the 
number of pixels across the telescope: 



tan© = 



D 

NH 



(6) 



and for a pixel number as assumed in Eq.[3] the angular size of 
one image element is: 



= arctan 



(7) 
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For a 10 km detection (at 532 nm) this amounts to « 0.8", and 
for a 30 km detection the angle seen is 0.4". These are within 
the size range of a seeing limited PSF. Additionally, global at- 
mospheric tip-tilt is completely excluded from the detection 
due to Fermat's principle. Therefore the image will be smeared 
out less than the seeing limit. In closed loop operation with 
high order AO, the situation will improve since the image will 
be sharpened. Larger scattering heights (e.g. sodium layer ex- 
citation) will result in a smaller maximum number of pixels; 
but due to the lack of cone effect there is no real advantage in 
going to greater heights. 

If the pixel number is chosen in the range as stated here, 
there is a perfect coincidence with the turbulent structure of the 
atmosphere: at short laser wavelengths more fringes are cre- 
ated, thus more detection pixels are required, and the turbu- 
lence is sampled on smaller scales. The ^128-pixel limit for a 
30 km, 532 nm detection would divide the telescope into sub- 
apertures of 6.25 cm squared, which is just the number required 
to sample the atmospheric turbulence on scales of Fried's co- 
herence length in the visible. Of course most projects will not 
need the high sampling that is required for short wavelength 
correction and it should be noted that the method outlined here 
can be used with fewer detection elements and a lower num- 
ber of fringes across the telescope. By changing the number of 
fringes, the sensitivity can be adapted to the wavelength of cor- 
rection and could even be chosen on-line to match the ambient 
atmospheric conditions. 

There is therefore a wide range of possibilities to arrange 
the gating height to fit the needs of an individual project: mul- 
tiple imaging over the propagation, using Rayleigh scattering, 
or sodium resonance excitation. 



occur is then only 0.0017 and decreases rapidly to 5 • 10~ 5 at a 
distance a chosen to be 0.8ro. 

Recalling Eq.[2the interferogram contains, besides the at- 
mospheric component, the phase term from the tilt O, of the 
laser wavefronts. If the standard PSI algorithms mentioned 
above are used, the <E>, component will always be present in 
the measurement, steadily increasing over the measured frame. 
This causes pixels on the detector to be close to the ambigu- 
ity of the tangent function, and with noise present in the mea- 
surement, 2n errors could easily occur. Removing those errors 
via phase unwrapping is computationally demanding. To avoid 
pixels on the wavefront sensor giving a measurement near the 
phase ambiguity and therefore being sensitive to noise, we pro- 
pose to use the phase retrieval algorithm outlined in the fol- 
lowing paragraph. The algorithm uses calibration frames to re- 
move the <b, component from the measurement. With two refer- 
ence frames and two on-sky measurements a solution is found 
from Eq.^2 Combining more phase shifts and reference inter- 
ferograms leads to an optimum phase measurement, explicitly 
shown for three steps in Eqs. 1121151 This treatment shifts the 
measurement for each pixel into the center of the unambiguous 
phase range and hence removes any need for phase unwrap- 
ping. 

If we assume that in advance of the on-sky measurement 
a set of interferograms is recorded that do not contain the $> a 
component, we get in total N s equations with i = {1,2, 3, ...N s ) 
for the N s different phaseshifts that are applied: 

/,<*, y, a) = I'(x, y) + I"(x, y) cos (<D fl (x, y) + <D,(x, y) + or,-) (9) 
I ci (x, y, a) = I'(x, y) + /"(*, y) cos (*,(*, y) + «,) (10) 



2.3. Phase retrieval algorithms 

There are several standard algorithms to calculate <I> fl from a 
set of phase shifted interferograms, which can be found in 
textbooks (e.g. Malacara 1992 1. With the four-step algorithm 
four interferograms U(x,y) are recorded with phaseshifts of 
or; = {0, | ,n, 3?} respectively. The solution to the four equa- 
tions contains only the measured intensities: 

h-h 



tan(0(x, y)) 



h-h 



(8) 



Comparable solutions can be found for three or more steps, 
e.g. the Hariharan or Carre algorithm. Common to all of them 
is that the resultant phase is expressed as the arctangent of the 
relation of measured intensities. Using modulo 2n correction, 
which is implemented in most computing languages as a two- 
argument tangent function, unambiguous results are obtained 
within a phase range of [—n, n]. Larger wavefront phases are 
usually calculated with the assumption of a steady function and 
a process called phase unwrapping. 

Since the method proposed in this paper only measures the 
phase difference between the sheared beams, phase unwrap- 
ping due to the sky component <& fl will not be necessary if 
the shear distance a is chosen to be of order ro. Assuming, as 
a back of the envelope calculation, a Gaussian distribution of 
gradients, the probability for phase differences larger than n to 



Id(x, y, at) are reference frames which are recorded with the 
same wavefront tilt, and the same phase shifts a, for which 
h(x, y, af) will be recorded on sky. The reference could either 
be taken by averaging lots of images from the sky, or a camera 
or special scattering screen could be inserted in the AO system 
conjugated to the height at which on-sky images will later be 
taken. 

The combination of Eg s . l9"l and 1 1 01 leads to an interesting 
solution for every consecutive i and i + 1 . By using trigonomet- 
rical identities the following equation can be found for each 
similar combination of 7y + j : 



O fl h - h\ +h- hi 4 ,ai-a 2 . 

tan( — ) = tan( ) 

2 h+hi-h-Ia 2 



(11) 



N s laser shots will give N s solutions for <£„ that are inde- 
pendent of the O, component and contain the difference of the 
two phase shifts that have been applied. It is noteworthy that 
<t> a is measured at half angle, which extends the phase range 
that can be directly detected to [—n, n], cancelling the need for 
a modulo 27r correction. In addition, there is no need to fix the 
phase shifts at a particular value, since they could be measured 
separately for the individual shots and just be included into the 
calculation. Explicitly for three steps the three solutions are: 

. ,<&a. h-Id+h-Ic2, .ai-a 2 I n l 

tan(— )i = — — — tan(— - — ) := — (12) 

I l\ + l c \ —12 — lc2 I Idl 
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tan(— ) 2 = — — — tan( — 

^ 12+ id - I-}- J C 3 



-as3 __ Ija 
2 hi 



(13) 



. ,^a, h-Ia+h-I c i or 3 

tan( — h = tan( — 

2 h+I c3 -h-I cl 



■a\. 



ha 

Id3 



) := -F (14) 



These equations show that, in principle, two laser shots and 
proper calibration frames could be enough to evaluate the phase 
difference. Due to the error behavior of the formulae above 
this is not the best choice. Each time the denominator becomes 
small, large errors occur in the calculation. If a weighted sum of 
the three (N s ) equations is used, a solution with minimal error 
behavior can be found, with the power k being an even number: 

1 

" l d\ 



tan(-^) = - 
2 /* 



d\ ^ 1 d2 ^ 1 d3 



hah-, + hih ' 



'd3 



(15) 



The algorithms will perform best if the a steps equal 2n/N s 
evenly distributed around a circle. The same principle can be 
applied for any number of a steps, and the detection becomes 
more stable against phase shifting errors. Due to the limited 
time in which the laser beams have to be launched, three or 
four seems to be a good choice. 

2.4. Detection errors 

The number of photons that is necessary per pixel to achieve 
a good measurement of d> fl is a crucial value, since it deter- 
mines the total laser power required. Each of the above men- 
tioned algorithms shows a distinct error behavior that depends 
on the choice of the phase shift a,, and for some algorithms on 
the actual phase value under test. In general it has been shown 
(Bruning 1992 1 that intensity fluctuations cause the standard 
deviation in the measured phase to be: 



1 



(16) 



where N s is the number of phase steps and S the intensity 
signal-to-noise ratio. In the case of pure photon noise the stan- 
dard deviation goes as (Wvant fT975> : 



cr 



I 

VP 



(17) 



where P is the total number of photons detected. Fig.[3]shows 
the variance of the phase gradients calculated numerically with 
the calibrated three- and four-step algorithms. In this figure the 
photon noise, and a contribution from readout noise and sam- 
pling of the interferogram, are included. Depending on the ac- 
tual scheme of detection and calibration, either algorithm could 
be preferable. With a total of 100 photons arriving within the 
three or four laser shots on one detector element, the cr 2 of the 
on-axis phase gradient can be as low as * 0.02rad 2 . 

Other error sources that might be special for this method 
are shot-to-shot variation of laser power and intensity distribu- 
tion. If this turns out to be a serious issue the actual intensity 
profile of the launched beam could be measured with a camera 
online and taken into account. Incorrect phase shifts are an- 
other possible error source that should be taken into account 
for a real system. While the limited coherence time of the at- 
mosphere calls for an algorithm with a low number of phase 
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Fig. 3. Error behavior of PSI algorithms. The top panel shows 
how the best achievable a 2 of the gradient measurement be- 
haves over the phase range that can be detected. For this plot 
100 photons per laser shot, a detector noise of two electrons and 
a sampling of n/4 has been assumed. Near the discontinuities 
of the tangent 2n, errors occur that appear as noise in this plot. 
While the noise in the four step algorithm (solid lines) does not 
depend on the actual phase, the three steps (dashed lines) shows 
a dependence on phase position. In the middle panel the depen- 
dence on photon number of the average variance is shown, in- 
cluding again a detector noise of 2 electrons and a similar sam- 
pling as before. In the bottom panel, the increase of cr 2 due to 
the sampling A from Eq.|3]is drawn, again assuming 100 pho- 
tons, but no detector noise. Up to a pixel size that spans « n/2, 
the noise in the measurement is only marginally increased. 



steps, the stability against phase shifting errors increases with 
the number of steps taken. Thus for a low number of steps an 
accurate phase shifting control should be present. With a fast 
electro optical phase shifter built into the system, this would be 
no major difficulty. 
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image polarization 
rotator rotator 

Fig. 4. Possible scheme for polarization combining in the two- 
direction gradient measurement. The linearly polarized laser beam is 
split in the interferometer into four overlapping beams. While two of 
these propagate directly to the PBS, the tilt and the polarization plane 
of the other two beams are rotated by 90° in the image and polariza- 
tion rotator. Combining all four beams is then done in the polarization 
beam splitter. On the way back the planes can be disentangled in front 
of the wavefront camera with a Wollaston prism. 



2.5. Gradients in two coordinate directions 

The sheared two beam detection, as mentioned above, only 
measures the gradient of the wavefront distortions along the 
direction of the applied tilt. To retrieve two-dimensional infor- 
mation additional means have to be employed. Apart from the 
possibility of disentangling x- and y-gradients via temporally 
independent measurements, two simultaneous methods will be 
described here: the use of orthogonal polarization states, and 
three beam interference. The first method is straightforward as 
shown in Fig. |4] one pair of tilted beams is launched linearly 
polarized. The second pair of beams has its tilt and the polar- 
ization state arranged orthogonally to the first pair. In a polar- 
ization beam splitter all the beams are combined and sent to 
the sky. Since the back scattering process is fully polarization 
preserving, the image of the scattered light on the wavefront 
sensor can again be separated with a polarizing beam split- 
ter into the desired x- and y-components. The calculation of 
the wavefront gradients can then follow fully the process de- 
scribed above. While the laser launching system in this scheme 
is quite simple and can be made fully polarization preserving, 
the path through the subsequent adaptive optics and image de- 
rotators of the main telescope may cause some trouble due to 
unwanted mixing of the two polarization states. A careful de- 
sign and eventually an active polarization control would be re- 
quired. 

Three or multiple beam interference can also be used to de- 
rive information along two directions. If three beams are coher- 
ently overlaid, two tilt directions can be applied. If the central 
beam is considered as number one, the second beam would be 
sheared and tilted along an x-direction and the third one along 
y. The resultant interference pattern is described by an equa- 
tion containing six unknowns, and therefore a total of six phase 
shifts will be required to solve for <S) ax and <£„,,. Depending on 
the speed of the pulsed laser and the detection camera this may 
be an option, but will take longer than the method described 



above. An alternative is the use of Fourier transform meth- 
ods to derive <S) ax and <t> flv . Methods exist to calculate gradients 
from multiple beam shearing interferograms with Fourier trans- 
form techniques (e.g. Primot[l993}. Shearing the three beams 
in the x- and y-direction and using phase shifts gives easy ac- 
cess to Fourier methods. In the three beam interference case, 
calculating the differences of two measurements — one with 
a x = a y = and one with a x = a y — n — will result in a 
sum of decoupled cosine functions in x and y. This can be fil- 
tered in Fourier space to disentangle the coordinate directions. 
A detailed analysis of this method is left for a forthcoming pub- 
lication. 

2.6. Sensitivity to turbulence at different heights 

The tilt of the coherent waves causes neighboring locations in 
one layer of the turbulence to be probed. If we consider first 
one probe of the phase, some basic relations of the sensing 
geometry can be found. At ground level, one phase will 'see' 
the atmospheric distortion at the point (x, y, zo), while the other 
traverses the point (x + a,y,zo). While traveling upwards the 
distance between the two paths decreases linearly, approach- 
ing zero at the height of detection. The original fringe pattern, 
formed due to the tilt, is therefore modified by the integral of 
the optical path length from zero to the detection height H. This 
can be regarded as a shear of the phase distortions by the dis- 
tance a at ground level, and a shear s at height z. 



s = a(H-z)/H 



(18) 



The measured phase difference <l> a is therefore most sensitive to 
the ground level part of the turbulence, since at higher altitudes 
nearly the same region of turbulence is probed. With differen- 
tial atmospheric phase distortions (f> atm super-imposed onto an 
on-axis light path, the measured signal is connected with: 



4>« = f < 

Jo 



„(x,y,z) - </> a tm(x + s,y,z)dz 



(19) 



For a subaperture, with minor high-order distortions, <t> a is di- 
rectly proportional to the average local wavefront gradient <f>' atm 
multiplied by the shear s: 



Jo 



• sdz 



(20) 



Due to the insensitivity to turbulent layers near the detection 
height, low altitude scattering will retrieve mainly the near- 
ground component. For example, if scattering at 20 km height 
is used, only half of the turbulence strength at 10 km can be 
seen. For sodium excitation at 90 km nearly all the structure of 
the turbulence is contained in the measurement. But even so, 
the remaining 10% of the distortions from a 10 km height will 
prevent diffraction limited correction from being achieved in 
the visible. 

It should be noted that in comparison with a 'standard' laser 
guide star, the phase disturbance of turbulence at higher altitude 
is contained within the measurement, albeit at smaller ampli- 
tude. The usual cone under a laser guide star totally misses the 
structure outside the cone, and does a wrong estimation of the 
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higher components, due to the projection geometry onto the 
pupil. 

To compare the possible performance of this scheme with 
a 'standard' laser guide star, we adopt as a first estimate the 
statistical treatment of the phase error, caused by focal aniso- 
planatism, as derived by Parenti and Sasiela [T994l 



^ = 0.5Z)f(|)^ec^- 0.904^] 



(21) 



for the part of the turbulence which is not measured outside the 
cone, and: 



^L, e = 0.057D!(^) 2 sec^ lH 



(22) 



for the part that is not sensed above the guide star. Here fi„ is 
the usual n-order moment of the turbulence structure function: 



(h)h"dh 



(23) 



For PSI sensing with tilted wavefronts, the weakly sensed — 
and therefore weakly corrected — part of the turbulence is pro- 
portional to the shear a - s(z). This assumption will hold as 
long as the distance a is not much larger than the atmospheric 
turbulence cells. The variance for the LPSI method due to the 
weakly sensed parts is therefore proportional to the first order 
ofyU„: 



.27T 



1 



= 0.057D3(— Y sec£-/i 



H 



Ho 



(24) 



The remaining phase error is always the combination of the 
unsensed part of the turbulence above the guide star, and a con- 
tribution of the insufficient measurement up to H: for single 
guide stars the part outside the cone, and for LPSI the under- 
estimation of the upper layers. For this estimation a simple 
Hufnagel- Valley model of the C\ distribution is assumed, as 
well as a power spectrum as present in Kolmogorov turbulence. 
The used Cl model is: 



Cl= 2.2 • 10" 23 /z 1 V' + 1.10 • 10 



+ 1.7- 10 



-14 - — 

1 ^ ,i n i 



(25) 



In Fig. |3 the variance of both possibilities is shown as a 
function of guide star (i.e. gating) height for an 8-m telescope. 
The linear dependence of the weakly sensed areas causes the 
LPSI method to perform better, especially in the mid-range 10- 
30 km detection heights. In total Strehl ratio (computed with 
the Marechal approximation), the difference for a telescope at 
this size is not really significant. For a ground layer adaptive 
optics system the linear decrease of sensitivity with height in 
LPSI could be a real benefit. Compared with single guide star 
systems, turbulence is never sensed wrongly or omitted outside 
a cone. In the following section we will show arrangements that 
allow for a full turbulence recovery. 

3. Multiple height measurements 

3.1. Gating Rayleigh scattered light at two heights 

By extending the proposed scheme to multiple gating heights, 
more information can be gained on the turbulence volume 




100.0 



Fig. 5. Variance of the phase against guide star height. Values 
are calculated for a 8.2-m telescope, with the C\ model from 
Eq. [25] assumed. Solid lines: single 'standard' guide star. 
Dashed lines: variance of the phase that is corrected with a sin- 
gle layer LPSI detection. The linear insensitivity to turbulence 
near the detection height shows an advantage over the quadratic 
dependence of the single guide star, especially if the scattering 
is at medium altitude. 



Fig. 6. Arrangement of the tilted waves when multiple scatter- 
ing from the atmosphere is detected. The camera is gated such 
that images from at least two different heights are retrieved. By 
choosing the right pixels, and subtracting them, all turbulence 
from below the first height can be recovered. The sampled tur- 
bulence volume is indicated with dashed lines. 



above the telescope. If we assume a scheme as shown in Fig.|6j 
the wavefronts are tilted such that a 'crossing' takes place at 
the first detection height. The resultant phase differences are 
detected at two heights at the points 01,2,^1,2, ••• correspond- 
ing to the pixel location of the detector. Taking now the two 
measurements of the phase difference at these heights and sub- 
tracting the shifted frames from each other, will result in a full 
on-axis measurement up to height H\ as can be shown easily: 
the measurement AO;,, at point b\ is the difference of the col- 
lected phase changes along the path aq — > b\ and bo — * b\\ 



4> 



(26) 
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Fig. 7. Solid lines: single 'standard' guide star. Dashed lines: 
Variance of the phase that is corrected with a double layer LPSI 
detection in dependence of scattering height h of the upper de- 
tection. The lower scattering plane is always chosen to be at 
half height of the upper one. Values are calculated as before for 
a 8.2-m telescope, with the C\ model from Ea.l25lassumed. The 
linear insensitivity to turbulence near the detection height now 
only shows an effect on the part of turbulence which is above 
h/2. At a detection height above 25 km the phase variance in all 
bands could be reduced below 0.2 rad 2 . 

and the measured phase difference at the point c 2 is: 



AO C2 = O fl0 , C2 - <D Co , C2 



= <D fl0iil + <5 ftljC2 - <D Co , Cl - <D CliC2 (27) 
Calculating the difference of the shifted frames results in: 



A<D C2 - A<5 6l = ® boM 



(28) 



The term <bb u c 2 - ^c u c 2 describes the path from Hi to H 2 . This 
measurement will show the decreasing sensitivity with increas- 
ing height to turbulent structures. The first term Q>b a M ~ ®c ,c, 1S 
the description of the full on-axis gradient in the height range 
from ground level to the first detection. The sampled path in 
the atmosphere is indicated in Fig.[6]with dashed lines. Similar 
relations can be derived for any use of two consecutive detec- 
tions. 

If the measured gradients could be fully turned into a cor- 
rection of the disturbed wavefronts, a system performance as 
shown in Fig.[7]would result. For this calculation the C 2 n model 
above has been used again. The detection height of the lower 
gating was assumed to be half of the upper height. Comparing 
again with a 'standard' single guide star, the system described 
here would give a much better performance, even at low scat- 
tering altitudes. If the upper detection is above «25km, the 
on-axis wavefront is sampled completely, so that even in the 
visible the remaining on-axis <x 2 is below O.lrad 2 . 

There are other possibilities to arrange gating heights and 
more detections. With several measurements the turbulent 
structure could be resolved in small height slices, leading to 
the possibility of assigning conjugated DMs to each of them; 
which despite the extra effort, is worth the gain. 

There may be some practical problems with the schemat- 
ics described above. Since we rely on the fact that the subtrac- 



tion of the two detections takes place in neighboring pixels, 
both should be made with the same number of pixels across 
the telescope aperture. Since scattering from the upper detec- 
tion results in a much lower photon flux and also the turbulent 
structure at larger height showing a larger coherence length, a 
sampling with larger subapertures would be desirable. Scaling 
low resolution upper layer data onto the highly resolved low 
altitude scattering might be possible, but would probably lower 
the performance of the system. 

As shown earlier, the sampling limit at a 30 km detection is 
around 128 x 128 pixels for a proper measurement at an 8-m 
telescope. If the same number is used for the lower detection, 
the subaperture size is x6 cm. This would be surely sufficient 
to correct in visible wavelengths and to gain high strehl in the 
near infrared. The only limit remaining is then the number of 
photons that can be produced in the scattering process. This 
number will be discussed later on. 

3.2. Multiconjugate correction with two beams 

The possibility shown above allows a full recovery of the tur- 
bulent structure of the atmosphere. If deformable mirrors with 
a large number of actuators are available, a correction can lead 
to high strehl ratios even in the visible. Nevertheless, the mea- 
surement performed with the above system is completely on- 
axis and will thus lead to the usual effects of angular aniso- 
planatism. In the visible the isoplanatic patch is rather small 
so only a small field of view can be corrected. By using multi- 
ple deformable mirrors conjugated to different heights, the field 
can be extended. 

The LPSI method is almost perfectly suited for this appli- 
cation. In Fig.[8]a possible system is shown. A first laser beam 
(consisting of overlaid coherent wavefronts) is directed over the 
first deformable mirror to the telescope, expanded and launched 
over the full aperture as described previously. This beam is de- 
tected at lower altitude and steers the first deformable mirror 
in a closed loop scheme. A second laser beam is then injected 
before that, hits the second DM and then joins the beam path of 
the first laser at a beam splitter. To enable the bams to be prop- 
erly separated, the two lasers could operate at different colors. 
The second laser is detected at higher altitude, and DM num- 
ber two is controlled, independently from loop number one, by 
this signal. This leads to a correction scheme that is compu- 
tationally not demanding, somehow similar to the 'layer ori- 
ented' concept proposed by Ragazzoni 2000 Treating the pro- 
cess mathematically, the lasers sample differential phases <f>' atm 
according to the shear s: 



®a2 



= ( H 
JO 

= f 

JO 



a(H v 



^dz 
^dz 



(29) 



Since beam number two is sent via the first DM to the sky, the 
actual phase which is sensed with it is the difference of <S> a2 and 
<D a i, and the laser wavefront leaves the telescope already cor- 
rected for <D fl i. This leads to following expression when setting 
additionally H 2 = M ■ H\ : 

Ma -a C m , a C H1 , 

<5>ca-<S>ai = MH ^ J q 4>' alm zdz+— <p atm (H 2 -z)dz(30) 
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Fig. 8. Schematic arrangement of a multiconjugate correction 
with lasers and LPSI. Laser 1 is directed over the first de- 
formable mirror, and the backscatter detected at lower altitude. 
Laser 2, operating at a different color, is directed first over the 
second deformable mirror, which is conjugated to a high al- 
titude turbulence. After passing a dichroic splitter it joins the 
beam path of laser 1 and the backscatter is detected from a 
higher altitude H 2 - 



How strongly this signal is applied to the DM is a function of 
the conjugated position, where it is located. The total correction 
over both DMs is therefore written as: 



(31) 



with Q being the wavefront shear at the location of DM2. If Q 
is chosen to be: 



Q 



M 



M- 1 



(32) 



all z-dependent terms in above equation cancel up to H t . The 
remaining terms are: 



= a Ktmte + Qjr ¥at 

Jo "2 Jh, 



SMi-z)dz 



(33) 



The first term describes the complete sampling of all gradients 
up to height H\ . The second term shows the linear decrease in 
sensitivity from H\ to Hi- The statistical treatment leads to a 
picture which is similar to the one shown in Fig. Sensing 
again at 10 and 20 km the turbulence is sampled nearly com- 
pletely — leading to possible high strehl ratios in all wave- 
length bands. It is noteworthy that in this scheme there is 
no computational subtraction of the two detection heights in- 
volved. The difference in the above equations results from the 
optical setup. Therefore it should be easily possible to sample 
the second laser with fewer subapertures. And since the control 



loops of the two DMs can be kept separate, sampling at differ- 
ent speeds of both lasers is easily possible. This would greatly 
relax the requirements for power of the second laser. 

It should be noted that one is completely free to choose 
the type of scattering involved. Both upper and lower altitude 
detection could be Rayleigh-type, or the upper scatter process 
could take place in the earth's sodium layer. 

Expanding the scheme outlined above allows one to gain 
additional information: when the lasers are launched slightly 
defocused, higher altitude layers are sampled outside the pupil 
of the telescope. Applying this information to larger de- 
formable mirrors would increase the off-axis correction. The 
phase retrieval will not suffer from such an expanded beam, 
since only relative phase differences are measured against a cal- 
ibration. The achievable corrected field depends on the struc- 
ture of the turbulence in the atmosphere, and how well the sec- 
ond DM is conjugated to a strong high altitude layer. A detailed 
analysis of this topic is out of the scope of this paper and will 
be given in a forthcoming publication. 

4. Photon numbers and laser requirements 

The number of photons required to achieve a good phase mea- 
surement is no different to the usual calculation for laser guide 
stars. The sensing geometry is just inverted: with a standard 
guide star the telescope pupil is divided into sub-apertures; with 
the LPSI setup, the full aperture of the telescope is used for the 
backwards imaging, but the pattern is divided by the number of 
image elements of the camera. In the following discussion, the 
photon number is counted to be the average over the fringes. 

For Rayleigh scattering the photon number N p h is calcu- 
lated from the standard theory by (see for instance van de Hulst 
1981): 



Nph = T] 



ED 2 p(H)^AH 



N 2 yH 2 



(34) 



E being the energy of the pulse, y the energy of one photon, 
p(H) the number density of molecules, H the scattering alti- 
tude, D the diameter of the telescope, N the linear number of 
pixel across the aperture and AH the length over which the scat- 
tering is sampled. r\ denotes the overall efficiency of a round 
trip from laser to receiver. is the scattering coefficient per 
unit solid angle £2, in the direction a with the angle (3 between 
scattering plane and polarization direction: 



dcr(a) 
dQ. 



(no- iy 



(cos 2 (yS) cos 2 (a) ■ 



• sin 2 (/3)) 



(35) 



no denotes the refractive index and A^o the number of particles at 
sea level. Using the usual density models of the atmosphere N p f, 
can be computed. For the plots in Fig|5]the USSA-1962 model 
from McCartney 1976 has been used. It should be emphasized 
that the number of photons incident on one detector element 
for a given laser power does not depend on telescope size but 
only on the chosen size of the subaperure D/N. 

With 'normal' gating of the wavefront camera, the scatter- 
ing altitude AH has to be kept very small to avoid a smearing 
of the image while the laser pulse is traveling this distance. The 
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Fig. 9. Photons collected by a single subaperture during one 
laser pulse with Rayleigh scattering in the atmosphere. The 
subaperture size is assumed to be 15 cm and the overall effi- 
ciency is taken to be 25%. Left panel: the situation for a gating 
AH of 90 m at 10 km and 280 m at 30 km are shown. Right 
panel: photons collected if a refocussing device were used. 
Refocussing would take place over 20% of the gating height. 
The laser wavelength is assumed to be 527 nm. 



angular size 6 of AH as seen from the outermost areas of the 
detector is: 



tan £ = 



DAH 



H 2 - HAH + D 2 



(36) 



which is equivalent to the focal depth of the camera. If we allow 
this angle to extend over half a pixel and combine Eq. [36] with 
Eq.[6]we find: 



AH = 



1 H 2 + ND 2 

2 H(l+N) 



(37) 



To avoid intensity being recorded in the wrong subapertures 
the possible gating height for an 8-m telescope at 10 km is only 
~90 m, increasing to 280 m at 30 km. This is a fairly small vol- 
ume out of which scattered light can be collected, and decreases 
even more for a larger telescope. This situation is reflected in 
the upper plot of Fig. |9] For a measurement with 15 cm sub- 
apertures, the photon flux from a 10 km detection height is just 
sufficient with a 20 mJ laser (20 W average power at 1 kHz rep- 
etition rate), when assuming an overall efficiency of 77=0.25 and 
a laser wavelength of 527 nm. A 30 km detection would require 
larger subapertures and very high pulse energy. 

The situation changes when using a device that is capable 
of keeping the image in focus during a part of the travel time. 
Proposed by Angel 2000 dynamic re-focusing would collect 



many more photons. Depending on the time during which the 
image could be stabilized, the power demands for the laser de- 
creases rapidly. This is shown in the left panel of Fig. [9] If the 
range over which photons are collected is increased to 20% of 
the scattering height a 20 mJ laser at 527 nm would be sufficient 
for a 30km detection with 15 cm subapertures. For the case of 
an extremely large telescope in the size range of 30 m or more, 
a refocussing device would be mandatory, since the focal depth 
of normal gating will decrease to very small values. 

Use of UV-lasers would be another possibility to increase 
the number of collected photons, due to the A dependence of 
the scattering coefficient. This has to be balanced against the 
quantum efficiency of available detectors, the optical efficiency 
of the system at UV-wavelengths and the available laser power 
at the desired repetition rate. 

Sodium layer excitation might be an option to collect pho- 
ton numbers comparable to a 30 km height detection more eas- 
ily. Due to the low repetition rate that is needed (^1 kHz, com- 
pared to usual 30-60 kHz for sodium guide star creation) sat- 
uration would become a serious issue in a small laser spot. 
Fortunately the area which is illuminated is very large. For an 
8-m aperture we calculate that saturation will not play a ma- 
jor role for a linewidth of 1 GHz with up to 20 W laser power. 
For extremely large telescopes the power density will decrease 
even further. 

Considerations concerning the beam quality of the laser are 
totally relaxed in the proposed system, compared with 'stan- 
dard' guide stars. For the usual sodium layer guide stars a beam 
quality of M 2 < 1.3 is usually mandatory. With the splitting 
and coherent overlap of the laser wavefronts and the calibrated 
LPSI algorithms proposed here, the initial wavefront distor- 
tions will cancel out completely and spot size is not an issue. 
In fact a typical top-hat beam profile of a powerful pulsed laser 
would be even better than a Gaussian profile from a TEMoo 
beam since the latter would under-illuminate the edges of the 
pupil. Lasers that fulfil these requirements are nowadays stan- 
dard products, available at reasonable cost. 

5. Extension to measure tip-tilt from the laser 

As already suggested by Ragazzoni 1995 the possibility exists 
to expand the area in which a suitable tip-tilt star can be found 
by using a small auxiliary telescope. It is worth mentioning that 
this ideally suited to the laser launching method proposed in 
this paper. From the laser pattern itself and the use of the main 
telescope for detection, the tip-tilt is perfectly invisible due to 
Fermat's principle. Therefore an additional reference star has 
to be used, as in any 'standard' LGS system. This puts con- 
straints onto the sky coverage of the system. Some areas which 
are of particular interest, such as deep fields, are selected due 
to the absence of bright stars. The use of the whole telescope as 
a laser projector, as proposed here, facilitates the use of an aux- 
iliary telescope for tip-tilt sensing. The principle is illustrated 
in Fig. [17]] With the small telescope capable of moving around 
over some area, a far off-axis star can be selected to be along 
the line of sight of the laser pattern projected onto the sky from 
the main telescope. At a height above all main turbulent layers 
of ~20-30km, the global movement of the scattered laser pat- 
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Fig. 10. Schematic showing how to measure tip-tilt above the 
main telescope by comparing the motion of the laser pattern as 
seen from a side telescope to that of an off-axis bright star. 

tern will represent the tilt component of the wavefront, since 
the laser beam has probed the full pupil size of the telescope. 
Subtracting the motion of the star from that of the laser patch, 
as seen from the side telescope, will result in a measure of the 
real tip-tilt along the science beam. For telescopes of smaller 
size the gain in sky coverage can be great; but for the case of an 
extremely large telescope the gain in sky coverage with the ex- 
ternal tip-tilt tracking might only be marginal, since the global 
wavefront slope can be measured with very faint natural stars. 

6. Conclusions 

We have proposed a method to measure and correct atmo- 
spheric distortions for astronomical telescopes using a laser to 
probe the turbulence. The use of laser wavefronts over the full 
aperture and the technique of phase shifting interferometry al- 
lows the distortions to be measured and corrected without suf- 
fering from the cone effect. The system outlined is therefore ap- 
plicable to any size of telescope including extremely large aper- 
tures. While the number of actuators on the deformable mir- 
ror and the number of detection pixels on the wavefront sensor 
have to be adopted accordingly, the laser system itself would 
not change with telescope size. We have shown that algorithms 
for the interferometric phase retrieval exist and allow for a cal- 
ibration of the measurement. The basic scheme outlined, in- 
volving a single detection of the scattered light, gives easy ac- 
cess to ground layer adaptive optics correction. By using of 
two deformable mirrors and two laser beams, a multiconjugate 
correction can be achieved. The necessary lasers are commer- 
cially available, and the power demands can be reduced with a 
dynamical refocusing device. The total system can be built rel- 
atively compactly, and the laser plus adaptive optics would plug 
into the telescope as one instrument. This removes the need for 
separate laser projectors, simplifying the implementation. The 
two examples given here of how to arrange the collection of 
the interferograms and the correction schemes, can easily be 
modified to fit the needs of an individual project. While the 
basic feasibility has been shown in this paper, further study is 
required in order to address detailed questions arising from the 
implementation or extension of such a system. Examples are 
field-of-view considerations in the multiconjugate case, studies 



for stray light suppression onto the science camera, and multi- 
color sensing for phase range extension. 
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